Dynamics of domain wall networks with junctions 
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We use a combination of analytic tools and an extensive set of the largest and most accurate three- 
dimensional field theory numerical simulations to study the dynamics of domain wall networks with 
junctions. We build upon our previous work and consider a class of models which, in the limit 
of large number N of coupled scalar fields, approaches the so-called 'ideal' model (in terms of its 
potential to lead to network frustration). We consider values of N between N = 2 and N — 20, 
and a range of cosmological epochs, and we also compare this class of models with other toy models 
used in the past. In all cases we find compelling evidence for a gradual approach to scaling, strongly 
supporting our no-frustration conjecture. We also discuss the various possible types of junctions 
(including cases where there is a hierarchy of them) and their roles in the dynamics of the network. 
Finally, we revise the Zel'dovich bound and provide an updated cosmological bound on the energy 
scale of this type of defect network: it must be lower than 10 keV. 
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I. INTRODUCTION 

As the early universe expanded and cooled down, it is 
believed to have gone through a series of phase transi- 
tions, at which networks of topological defects must nec- 
essarily have formed [l|. The type of defect that forms 
and its specific properties depend on the particular de- 
tails of each symmetry breaking, and hence there is a 
wide range of possibilities, which will lead to correspond- 
ingly very different cosmological consequences [2]. Do- 
main walls are known to be pathological except if they 
are very light 3], but it has been claimed [J] that if a 
domain wall network is frozen in co-moving coordinates 
(or 'frustrates', as is often colloquially put) then it can 
naturally explain the observational evidence that points 
to a recent acceleration of the universe. 

It is clear that this scenario is subject to a number of 
observational constraints. For example the cosmic mi- 
crowave background (CMB) data [5| severely constrains 
the characteristic scale of the network, L, which needs 
to be tiny in order not to give rise to exceedingly large 
CMB fluctuations. Also, recalling that the equation of 
state of a domain wall gas is given by 



the walls and that we require that 
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with v being the root-mean squared (RMS) velocity of 



'Electronic address: ppavelin@fc.up.pt 



^Electronic address: Carlos.Martins@astro.up.pt 
t Electronic address: jmenezes@fc.up.pt 
^Electronic address: rms@fisica.ufpb.br 
^Electronic address: jeolivei@fc.up.pt 



W<-\\\ 






(2) 



in order to accelerate the universe at the present time 
it is clear that only a non-relativistic value of the veloc- 
ity (in fact, basically v ~ 0) would have any chance of 
working. The simplest domain wall models are known to 
reach a scaling regime (until they dominate the energy 
density of the universe), as first pointed out in 6] and 
recently studied in detail in [Tj, |8|, [9( , and hence are ob- 
viously unable to satisfy these constraints. Nevertheless, 
it was thought that more complicated models, notably 
those having junctions, would eventually frustrate and 
therefore might conceivably be able to satisfy them. 

In previous work [lfj . [ll| we have studied the dynam- 
ics of domain wall networks with junctions, and investi- 
gated in detail energy, geometrical and topological con- 
straints on their properties. This led us to develop an 
ideal class of models which includes, in the large N limit, 
what we called the ideal model (that is, the best possi- 
ble candidate for frustration). A series of analytic and 
numerical arguments then led us to a no frustration con- 
jecture: even though one can build (purely by hand, as 
was done in [TlJ, [t2|, LL3|) special lattices that would be 
locally stable against small perturbations, no such con- 
figurations are expected to ever emerge from any realistic 
cosmological phase transition. (A much simpler analysis 
of the evolution of non-interacting and entangled cos- 
mic string networks using a velocity-dependent one-scale 
model had already led to a similar result, at least in four 
space-time dimensions [14;].) Our high- resolution numer- 
ical simulations of various realizations of the ideal class 
and other models showed clear evidence of a gradual ap- 
proach to scaling, which was subsequently confirmed in 



[l5| (though note that the latter only reports on numer- 
ical simulations in Minkowski space thus neglecting the 
role of the expansion) . 

Still, a fair criticism of our earlier work [10, [ll( would 
be that we only considered numerical simulations in two 
spatial dimensions, and the behavior could well be differ- 
ent in the three spatial dimensions of our universe. More 
recently [161 ] we endeavored to eliminate this shortcom- 
ing and presented results of a series of massively par- 
allel domain wall network numerical simulations of the 
'ideal' class of models in three spatial dimensions, which 
confirmed our earlier work, providing conclusive evidence 
for a gradual approach to scaling and hence strongly sup- 
porting our no frustration conjecture. The present article 
is the fourth in this series of papers, and its aims are to 
present the results of [ill more thoroughly, to extend 
them in a number of different ways, and to discuss in 
some detail the corresponding cosmological implications. 

In Sect. [IT] we introduce various approaches to the 
study of domain walls — the microscopic description, the 
analytic modeling and the building of phenomenological 
models — that we will be using in the rest of the paper. 
In particular we highlight some relevant scaling solutions 
for domain wall networks and describe the ideal class of 
models which will be our primary focus. In Sect. IIIII we 
discuss several technical details of our massively paral- 
lel field theory numerical code (which uses an improved 
version of the algorithm of Press, Ryden and Spergel 6], 
henceforth referred to as the PRS algorithm). In par- 
ticular we describe in some detail how we identify the 
domain walls and how we measure their velocities. Sect. 
IIVI contains our main numerical results: we describe the 
outcome of our various sets of three-dimensional high- 
resolution simulations of the ideal class of models. We 
also discuss the relevance of the key findings and contrast 
the results with those of other models. This then leads 
us to discuss, in Sect. [V] the analytic modeling of our 
results, including the cosmologically important asymp- 
totic limit of the ideal class of models. In Sect. IVII we 
present an analysis of the possible hierarchies of junc- 
tions in the ideal and other models, as well as a brief 
discussion of the role of the junctions on the dynamics of 
the networks. Finally Sect. IVIII has our conclusions, as 
well as some comments on the implications of our results 
for cosmological scenarios involving domain walls and a 
brief outline of future endeavors in this area. 



II. MODELING DOMAIN WALLS 

We shall start by describing some relevant micro- 
physical and averaged (macro-physical) tools for mod- 
eling domain walls and their evolution. While parts of 
this section are a review of previous work (see for exam- 
ple (2J), they will be useful not only as a means to set 
up some notation but also as a basis upon which we will 
build the discussion in subsequent sections. 



A. What is a domain wall? 

Let us start by a basic but important point: defin- 
ing what a cosmological domain wall is in a simple but 
relatively generic and physically meaningful way. A cos- 
mological domain wall is effectively a two-dimensional 
object, meaning that its thickness is much smaller than 
the curvature radius. This implies that locally the do- 
main wall is planar. Consider a local inertial reference 
frame in which the domain wall is instantaneously at rest 
and assume that the domain wall is perpendicular to the 
z direction. Another defining characteristic of a cosmo- 
logical domain wall is that its properties are independent 
of x and y, which means that both its thickness, 5, and 
its mass per unit area, a, do not change along the wall. 
The fact that nothing changes along the wall implies that 
the physical velocity of the domain wall is perpendicular 
to it (it is not possible to measure parallel velocities). 

These simple properties are sufficient to determine 
their energy- momentum tensor |17|. Let us start by 
considering a domain wall network which is frozen in 
co-moving coordinates in a Friedmann- Robertson- Walker 
(FRW) universe, and let us assume without loss of gen- 
erality that the characteristic length is small enough 
that the network can be considered homogeneous and 
isotropic on cosmological scales. If that is the case then 
the average domain wall density is 



P = 



oA 
V 



(3) 



where V is a large volume and A is the total domain wall 
area inside V. Using the fact that 



dp 
~dt 



+ (1 + w)Hp = . 



(4) 



we get that w = —2/3 for a frozen domain wall network. 
Let us now move from the domain wall network as a 
whole and try to compute the energy momentum tensor 
T Q/3 of the individual domain wall segments. Consider a 
local inertial frame in which the domain wall is instanta- 
neously at rest and assume that the wall is perpendicular 
to the z direction. In this case, we can show that 



T„ 



diag [p w (z), -p w {z),-p w {z), 0] 



(5) 



where p w is the wall energy-density and Tu are the stress 
components along the i direction. The fact that the wall 
properties do not change along it means that it must be 
impossible to measure velocities in any direction parallel 
to the wall plane. So T^ v must be the invariant with 
respect to Lorentz boosts in any direction along the wall. 
Consider the following boost in the x direction 

t'= r y(t + vx), x' = j(x + vt), y' = y , z' = z (6) 

1/2 

where 7 = (l — w 2 ) is the Lorentz factor and v is the 

wall velocity. Hence T^ v transforms as 



T"'"' = k£k v pT a(s 



(7) 



where 

Af=Af= 7 , a£=A*'= 7V , A^'=Af = l (8) 
and all other components vanish. By requiring that 



rpt'v' — jj>ty _|_ y V j<xv — j> ty 



and 



T x ' v ' = jT xy + ryvT ty = T xy 
we obtain T ty = T xy = 0. Similarly from 



rpt'z' _ ^ T tz + ry V x xz = T tz 



and 
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we have T* z = T xz = 0. On the other hand, 
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implies that T xx — —T a . From the fact that there is 
no preferred direction along the wall we have T xt — T yt , 
T xz = T y Z ^ and TXX = T y V ^ FmaUy from the en ergy- 

momentum conservation law we have, T^ a a — 0, and 
so T z \ z = 0. This is equivalent to T zz = since T zz 
is null, at z — > oo, for any domain wall solution. We 
then conclude that a planar domain wall solution has the 
energy-momentum tensor given by Eq. [5l Consequently, 
given T 00 all other components of the energy-momentum 
tensor of the domain wall are uniquely determined. Note 
that all of this is independent of the microscopic structure 
of the domain wall. 

This has an important consequences for domain wall 
dynamics. Given the above considerations, a domain wall 
is essentially characterized by its mass per unit area a 
(which is fixed) and its velocity. Therefore there is only 
one variable to evolve, v, and its evolution is fully deter- 
mined by energy-momentum conservation T a ^ .p = 0. In 
the particular case of a planar domain wall one has 
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Now consider the role of friction effects on the dy- 
namics of domain walls. Since we are mainly interested 
in showing that no realistic domain wall network will 
frustrate by the present epoch we need only consider 
the maximum frictional force consistent with energy- 
momentum conservation. For interaction with radiation 
one has 



dv 

u 

while for matter we have 

' dv s 
It 
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(15) 
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and in both of these we are neglecting relativistic correc- 
tions. The extra factor of v which appears in the interac- 
tions with matter gives us an extra reason not consider 
them in the analysis. In our previous work [lOj, |ll| we 
justified not considering interactions with matter on the 
grounds that they would adversely affect structure forma- 
tion. However, there could still be some fraction of the 
dark matter in the universe (say one part in a hundred) 
with which domain walls could interact without causing 
obvious observable effects. However, that interaction in 
the small v limit (the one we are interested in) would be 
much weaker than in the radiation case. This justifies, 
and indeed strengthens our conclusions. 



B. Generalized Lagrangians 

In the above we considered some properties of domain 
walls which are independent of the particular scalar field 
model we were considering. But are the above conclu- 
sions regarding domain wall dynamics valid more gen- 
erally? It turns out that a more general argument can 
indeed be given. 

Consider a single-field model 



S= dt d 3 xyf^ £((/), X) 



(17) 



Varying the action with respect 



where X = 1/2<9 P 

to 4> we get the following equation of motion 

1 
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(18) 



where Cx = §§? and C^ — ^§. The FRW metric has 
been assumed. Re-writing this in physical coordinates 
we have 

and for planar domain wall solutions (depending only on 
one spatial coordinate, z) this can be re-written 



9 (c 9 ^ + w d<j) d (c dct, \-r 



which has a solution 
this reduces to 



(20) 
with 9 = z. For a static wall, 



dz \ dz 



C d 



(21) 



Now consider a wall moving with speed v in the z direc- 
tion. Then 
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(22) 
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noting that in this case X - — h{^ ) 2 , one can write 
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3H(vj) 



= £$ (24) 



and given that 

have 



is a solution to Eqn. 1211 we must 



d(vy) 
dt 



3H(v-f) = 



(25) 



which again leads us to jv cx a -3 , generalizing the result 
in the previous sub-section. 

In models with multiple scalar fields there are other 
ways of putting together derivative terms which preserve 
Lorentz invariance, for example d^cjjd^Xi but a similar 
analysis would still lead to Eqn. [25] There are, however, 
ways to circumvent this constraint. The most obvious is 
to relax the FRW metric assumption, or in other words 
moving away from standard cosmology. One example 
would be to consider models with extra dimensions. 



C. Simple analytic modeling 

In [7[ we introduced a phenomenological one-scale 
model for the evolution of domain wall networks that 
has been shown to provide a good approximation to the 
evolution of two key network parameters: the character- 
istic scale of the network, L, and the RMS velocity of the 
domain walls, v. Their evolution equations are 



dL 



L 



— = HL + -v* 

dt Li 



cv , 



dt y } \L l d 



(26) 



(27) 



where H is the Hubble parameter, c is the energy loss effi- 
ciency, k is the curvature parameter and we have defined 
a damping length scale, 



h= 3H+ r f 



(28) 



which includes both the effects of Hubble damping and 
particle scattering. The characteristic scale of the net- 
work is defined as L = a / p where p is the average den- 
sity in domain walls and a is the wall mass per unit area. 
Note that if domain walls are an important contribution 
to the dark energy then p must be of the order of the 
critical density, p Cl at the present time. 

Just as in the case of cosmic strings, the generic at- 
tractor of the above equations in a cosmological model 
where the scale factor grows as a cx t a is a linear scaling 
solution (although in the case of the domain walls they 
will eventually dominate the energy density of the uni- 
verse), that is characterized by the following (constant) 
parameters 



^ (L\-_ k(k + c) 
t J 3a(l — a) 



(29) 



3a k + c 



(30) 



If we ignore the energy loss by the network (by making 
c = 0) it is easy to show that a linear scaling solution is 
still possible for a > 1/4. In the radiation era we obtain 



L = —=kt 
\/3 



while in the matter era we have 
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L 



hi 
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(31) 



(32) 



(33) 



(34) 



Notice that in both eras we have L ~ kt and relatively 
large velocities. If we require CMB temperature fluctua- 
tions generated by domain walls on scales of the order of 
Hubble radius to be smaller than 10 -5 then one would 
need (iff) 3 / 2 < 1CT 5 or equivalently L < 1 Mpc (our 
conservative estimate in [l(|)- However, CMB observa- 
tions imply that the fluctuations generated by the do- 
main walls have to be smaller than 10 -5 down to much 
smaller scales (say ~ i/ _1 /100). This means that cur- 
rent constraints on L are expected to be roughly 2 orders 
of magnitude stronger, 

L<10kpc (35) 

which implies a very small curvature parameter 

fc<l(T 6 . (36) 

This clearly shows that the simplest domain wall scenario 
without junctions (which necessarily has k ~ 1) is ruled 
out as a dark energy scenario, regardless of any other 
considerations. It is easy to show that allowing for a 
non-zero c leads to a larger L and consequently it does 
not help frustration [3, [l0| ■ On the other hand, as was 
pointed out above, including friction also does not help 
much, due to the limited amount of energy with which 
domain walls can interact conserving energy and momen- 
tum in. 



D. Towards the ideal model 

This failure of the simplest domain wall scenario led us 
to consider more complex scenarios with junctions and in 
|lfj | we set out to investigate in detail energy, geometri- 
cal and topological considerations that severely constrain 
the properties of domain wall networks. In particular, in 
the context of 2D domain wall networks, it was shown 



using local energy considerations that two-edge domains 
are always unstable and that three-, four- and five-edge 
domains will be unstable if only Y-type junctions occur 
in a given model. We have also demonstrated that in- 
creasing the average dimensionality of the junctions, (d), 
leads to a decrease of the average number of edges, (x), 
per domain (in particular if (d) > 6 then (x) < 3 and con- 
sequently no equilibrium configurations will ever form). 
Moreover, allowing for domain walls with different ten- 
sions contributes to increasing the instability (relative to 
models where all the walls have the same tension). This 
is because the walls with higher tension will tend to col- 
lapse, thereby increasing the dimensionality of the junc- 
tions which, in turn, will in general lead to the production 
of further unstable two edge domains. 

Another important aspect is related to the fact that 
the stability of a given domain depends on global consid- 
erations (those depending on the configuration of the sur- 
rounding domains) as well as local ones (those associated 
with the domain itself) and we expect the global ones to 
become more important as we increase the dimensionality 
of the junctions or consider specific domains with a large 
number of edges. In particular it is possible to show (in 
2D) that a domain with three edges only survives the lo- 
cal stability analysis by very little in the case where there 
are only X type junctions (the potential energy after the 
collapse would be at most about 10% larger than before). 
In this case, we expect that, in general, non-local effects 
will make a three edge domain unstable. As a result in a 
model with only X type junctions the only possible sta- 
ble configuration is the one in which all the domains have 
the same number of edges (namely four), which never 
occurs in the context of realistic domain wall network 
simulations. Note that we are assuming the junctions to 
be free throughout the paper. Otherwise their energy- 
momentum contribution could not be neglected, spoiling 
the dark energy properties associated with a static do- 
main wall network. We shall return to the role of the 
junctions later in this paper. 

By combining energy, geometrical and topological con- 
siderations, one is then naturally led to a class of mod- 
els with any number N scalar fields and correspondingly 
N + 1 vacua, with the property that all possible domain 
walls have equal tensions [ll|. The 'ideal' model, mean- 
ing the one with maximal a priori probability of reaching 
a frustrated state, is then obtained in the limit N — > oo. 
For large N the collapse of a single domain will only very 
rarely lead to the fusion of two of the surrounding do- 
mains. This is clearly a very desirable feature from the 
point of view of frustration (see [n|). Also, by requiring 
all domain walls to have equal tensions we avoid another 
potential source of instability. 

A specific realization of this class of models with N 
scalar fields and a scalar field potential with N +1 vacua 
was first given in [11| 



where pi . are the N + 1 coordinates of the vacua of the 
potential. We have chosen pi . to be the vertices of an 
(N + l)-dimensional regular polyhedron, and fixed the 
distance between the vacua to be equal to the parameter 
tq. Given that in this model all possible domain walls 
have equal tensions, only Y-type junctions will form. Al- 
though it is possible to have only X type junctions in a 
model with 4 minima, this is in general not the case for 
a larger number of minima if all possible vacuum config- 
urations arc allowed to exist. In particular this means 
that it is not possible to construct a generalization of the 
'ideal' model in which all junctions are of the X type. 

In [ll| we performed 2D field theory simulations of the 
'ideal' class of models with N = 4 and N = 7. These 
results are useful since they are much simpler (hence eas- 
ier to understand), containing many important features 
which are also relevant in higher dimensions. However, 
realistic 3D domain wall network simulations were re- 
quired in order to test our no frustration conjecture. The 
results of a set of 3D matter era simulations was first pre- 
sented in [lj|, confirming the above expectations. In the 
following sections we will describe the algorithm used to 
obtain these in somewhat more detail, and present and 
discuss the results of a much more extensive set of simu- 
lations, 



III. MASSIVELY PARALLEL SIMULATIONS 

We have performed high-resolution field theory nu- 
merical simulations of domain wall networks on the UK 
Computational Cosmology Consortium's COSMOS su- 
percomputer using a modified version of the PRS 
Earlier descriptions of our code can be found in 
and some more specific details are described below. The 
numerical code, in C language, was parallelized with 
OpenMP directives, and gradually optimized for the 
shared memory architecture of COSMOS. The code plus 
its auxiliary routines are about 5000 lines long. 

The results to be discussed below come mainly from 
simulations of 128 3 , 256 3 or 512 3 boxes, in various cosmo- 
logical epochs (matter domination, radiation domination 
and so on) for the ideal class of models with all values 
of N between 2 and 20. In addition, some realizations of 
other models were also performed (these particular mod- 
els will be described as they are introduced) . Unless oth- 
erwise is stated, we assume initial conditions where the 
scalar field at each point in the grid is associated with a 
randomly chosen minimum of the potential. 

The required memory for our three-dimensional simu- 
lations is approximately 



MEM = N * 4.5 * 8 * 



\ 1024 ) 



GB: 



(38) 
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((pi-Pi,) 2 , (37) 



as an example, for 512 3 boxes 90 GB are required for 20 
fields. An output box binary file can also be produced 
at specified time steps which can then be used to 



generate animations: an example of these is available at 
http : //www . damtp . cam . ac . uk/cosmos/viz/movies/evo 
(several others are available from the authors upon re- 
quest). As a benchmarking example, a 512 3 simulation 
with 3 fields (requiring about 14.5 Gb of memory) takes 
about 14 seconds per step on 16 processors and only 5 
seconds per step on 32 processors (as the memory ratio 
becomes favorable) without box output, and a complete 
run takes just 85 minutes. For larger runs the scalability 
is good if one keeps the memory smaller than 1 GB per 
processor. The example given above, a 512 3 box with 20 
scalar fields and box outputs at every time step, takes 
about 2 hours and 15 minutes on 128 processors. 




A. Identifying the walls 

We recall that the key idea behind the PRS algorithm 
Q is to modify the domain wall thickness in order to 
ensure a fixed comoving resolution, thus enabling much 
larger dynamic ranges to be simulated. We defined the 
domain wall as the region where 



V{4>) > aV m 



(39) 



Knax being the maximum of the potential and < a < 
1. The chosen threshold value a will thus define two 
properties of the network: the corresponding thickness of 
the static domain wall 6 and the volume fraction of the 
box with domain walls /. They can be used to estimate 
the comoving characteristic scale of the network 






(40) 



where / is the volume fraction of the box with domain 
walls (that is, satisfying the above condition), and 8 is 
the thickness of a static domain wall. 

We have verified that for 0.2 < a < 0.6 (see FigHJ, 
our results are almost independent of the threshold a, as 
long as the domain wall is sufficiently resolved. 

In order to check whether the domain wall network is 
evolving according to a linear scaling solution we define 
the scaling exponent A so that 



ex 77 



,i-A 



(41) 



where r\ is the conformal time. If A = we have linear 
scaling, whereas the network frustrates for A = 1. The 
scaling exponent can be calculated at two different values 
of the conformal time 771 and 772 by using the following 
relation M 



Hm,m) 



ln(fli/fl 2 ) 

In (771/772) 



(42) 



where R = n L~ l . In general, since the simulations are 
evolved until 77 becomes equal to one half of the box co- 
moving size (because beyond that point the effect of the 
finite box size will contaminate the results) the values 



FIG. 1: The relative comoving characteristic scale, L c /n, for 
various values of the parameter a defined in Eqn. 1391 From 
bottom to top these are respectively a — 0.6 (blue, dotted 
line), a = 0.4 (red, dashed line) and a = 0.2 (green, solid 
line). All of these were obtained for 128 3 simulations. 



of 771 and 772 are chosen in order to consider the second 
half of the dynamical range of the simulations. Note that 
the earlier part is not included to avoid contamination of 
the results by our particular choice of initial conditions, 
while beyond one half of the comoving size, the effect of 
the finite box size will contaminate the results. 



B. Measuring wall velocities 

We measure the domain wall velocities using an al- 
gorithm analogous to that described in [7| which sep- 
arates most of the radiated energy from the walls. We 
emphasize that this is mandatory for any reliable velocity 
measurement, since otherwise this radiation would con- 
taminate the estimates of the velocities — our previous 
work shows that the contamination can easily be at the 
twenty percent level. This algorithm therefore represents 
an important advantage over previous velocity estima- 
tions (such as the one in the original PRS algorithm [y]). 
To be specific, we estimate the velocities as 



(*V) 
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(43) 



where a dot represents a derivative with respect to con- 
formal time and 7 = (1 — v 2 ) 1 / 2 . 

Let us revisit our discussion in Sect. Ill Al showing that 
71* oc a~ 3 for any planar domain wall perpendicular to 
the zz direction and that our method for determining 
the RMS velocity of the domain walls will work inde- 
pendently of the particular scalar field potential we are 
considering. Let us write the solution for a moving wall 
as (f)(9) where 38/dt — jv and d6/dz = 7. Then we have 



i)t 



d9 1V dz 



cW 
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(44) 



and from the fact that (f)(8) (with 9 — z) is a static solu- 
tion we have 



d 2 cj> dV 
d9 2 ~ dd> 



we immediately obtain 

d(wy) 
dt 



■3H(wy) = 0, 



(45) 



(46) 



which has the solution jv oc a 3 as claimed. Note that 
by integrating 



d 2 (j> _ dV 
~W 2 ~~d4> 



(47) 



we obtain 



00 



2V. 



(48) 



and consequently it is straightforward to show that 

2 



v 2 ^ 



1 

2V V dt 



(49) 



independently of the model. Applying the PRS proce- 
dure is equivalent to making the substitution t — > rj, 
z — > q z in the above equations, where dz = a(t)dq z and 
dt = a(t)drj. 



IV. NUMERICAL SIMULATION RESULTS 

We are now in a position to present the bulk of our 
work: the most extensive, most accurate and largest (in 
terms of the conformal time dynamic range, which is the 
meaningful diagnostic) set of 3D field theory simulations 
of domain wall networks to date. We will start by de- 
scribing the simulations of the ideal class of models in 
the matter era (which is the relevant context for dis- 
cussing the possible role of domain wall networks as dark 
energy). Later in this section we also discuss simula- 
tions of the ideal class of models in other cosmological 
epochs, as well as some representative cases of simula- 
tions of other 'non-ideal' models: we will consider some 
examples from the phenomenological models of Kubotani 
[18j | and of Bazeia et al. [19| (to which we will henceforth 
refer as the BBL model). When appropriate we will also 
compare these results with those obtained for simulations 
in other dimensions [lCj, LUj as well as those obtained for 
the simplest model of domain walls 0, @, Q i which have 
no junctions. 



A. Ideal model, matter era 

In Fig. [5] we plot the scaling exponents, A, defined by 



V 



oc r\ 



(50) 



for the ideal class of models with all numbers of fields 
N between 2 and 20, in the matter-dominated era. We 
show results of simulations in 128 3 , 256 3 and 512 3 boxes, 
to illustrate the effect of the box size (and hence dynamic 
range) on the results — these effects are small but still 
identifiable, even with these very large boxes. The error 
bars represent the standard deviation in an ensemble of 
ten simulations. 

We see that A is slightly greater than zero which in- 
dicates that there are small departures a the scaling so- 
lution (keep in mind that exact linear scaling obviously 
corresponds to A = 0). In other words, the wall net- 
work is not equilibrating as fast as allowed by causality. 
Overall we find 



A 



id. mat 



0.11 ±0.05. 



(51) 



which we can compare with the value we obtained [7| for 
3D 512 3 matter era simulations of standard domain walls 



X, 



0.04 ±0.02. 



(52) 



The two values are in fact relatively similar; the larger ex- 
ponents in the case with junctions should be attributable 
(at least in part) to the longer time needed for the relax- 
ation to scaling. The fact that as we increase the box 
size, thus evolving the simulations for a longer dynamic 
range, A gets closer to zero is a clear indication that the 
networks are slowly approaching a scaling solution. On 
the other hand, frustration (which would correspond to 
A = 1) is clearly ruled out. 

In Fig. [3] we plot the asymptotic values of u* = jv and 
L c jr\ for the above 256 3 and 512 3 runs. (Again, the er- 
ror bars represent the standard deviation in an ensemble 
of 10 simulations.) As in the previous plots it is notice- 
able that the differences between the results obtained for 
boxes of different sizes are visible but quite small — in 
general the results are consistent with one another, given 
the error bars. The way the differences go is also what 
one expects. Mindful of these differences, we can still 
say that for many of our purposes the 256 3 or even 128 3 
boxes (each of which we can usually run in about an hour 
or less) will produce adequate results. 

As we increase N the asymptotic value of L c /rj de- 
creases, which is expected as we get closer to the (asymp- 
totic) ideal model. It is also significant that the differ- 
ences between the successive N results for L c /tj become 
increasingly smaller for large N which is a clear indica- 
tion that the results obtained for N = 20 are already 
reasonably close to the N — > oo results. On the other 
hand, we do not find any significant dependence the ve- 
locities with N. Indeed tipically wc find 



V ld ,mat = 0.36 ± 0.02 , 



(53) 



-1/2 



where we estimated v as u* = v(l — v")~'<" or equiva- 
lently v = v*(l + v 2 )^ 1 ' 2 . Again we can compare this 
with the value we obtained for analogous simulations of 
standard domain walls 



v s t,mat = 0.39 ± 0.02 : 



(54) 
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FIG. 2: The scaling exponents, A, for all TV's between 2 and 
20, for the 128 , 256 boxes (top panel, top and bottom set 
of lines respectively) and 512 3 boxes (bottom panel) in the 
matter-dominated era. The error bars represent the standard 
deviation in an ensemble of 10 simulations. Note that frus- 
tration would correspond to A = 1. 



thus the velocities are slightly smaller in the case of net- 
works with junctions though given the various numerical 
uncertainties the difference between the two is not very 
significant. However, the scaling density in the case with 
junctions is much larger than in the standard case — by 
a factor of about 3 for N = 2, and increasing with the 
number of fields (we will discuss this dependence later 
on). A more important comparison can be made with 
the value we would expect if the network had no energy 
losses. In this case (recall Eqn. [34j) we would expect 
v ~ 0.41, and the difference means that energy losses are 
still noteworthy, despite the existence of the junctions. 
This is ultimately the physical reason why the frustra- 
tion mechanism can never be realized. 

This difference between the behavior of the two aver- 
aged quantities is crucial — although changing the number 
of fields N will change the network's characteristic length 
(and therefore its density) , the fact that the velocities do 
not change is an indication that in the 'local' dynamics 
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FIG. 3: The asymptotic values of v* — *yv and L c /r] (dotted 
and solid lines respectively) for 256 3 (top panel) and 512 3 
(bottom panel) matter era simulations of the ideal class of 
models with N ranging from 2 to 20. The error bars represent 
the standard deviation in an ensemble of 10 simulations. 



of each individual wall until the collpase will not strongly 
depend on N although the probability that the collapse 
of a given domain will result in the fusion of adjacent 
domains will be higher for smaller N. Note that we are 
assuming that the junctions themselves are dynamically 
unimportant. 



B. Ideal model, other eras 

Even though simulations of the matter era (that is, 
with the scale factor evolving asooc t 2 / 3 ) are the most 
relevant from the point of view of cosmological scenar- 
ios involving domain walls, it is also interesting to study 
their evolution in other cosmological epochs. Results of 
analogous simulations of the ideal model for the radiation 
era (a oc i 1 / 2 ) are shown in Figs. 2] and [5] The scaling 
exponents and properties are defined and measured as 
before, and the simulation parameters are also similar, 
the only difference being that in this case we only have 
carried out 128 3 and 256 3 simulations. 
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FIG. 4: The scaling exponents, A, for all N's between 2 and 
20, for the 128 , 256 boxes (dotted and solid lines respec- 
tively) in the radiation-dominated era. The error bars repre- 
sent the standard deviation in an ensemble of 10 simulations. 
Note that frustration would correspond to A = 1. 
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FIG. 5: The asymptotic values of v* = fv and L c /rj for radia- 
tion era simulations of the ideal class of models with N rang- 
ing from 2 to 20 (top and bottom lines respectively). These 
come from simulations of 128 3 (top panel) and 256 3 (bottom 
panel) , and the error bars represent the standard deviation in 
an ensemble of 10 simulations. 



The results are qualitatively identical to the ones in the 
matter era, although as expected there are some quanti- 
tative differences in the scaling parameters. We now find 
the scaling exponent 



A„ 



0.10 ±0.05. 



(55) 



which is again comparable to the one obtained for anal- 
ogous simulations of standard domain walls 



Xst,rad = 0.04 ±0.02. 



(56) 



Likewise, we do not find any significant dependence the 
velocities with N. We now obtain 

Vid,rad = 0.45 ± 0.03 (57) 

whereas for standard domain walls in the radiation era 



V s t,rad = 0.48 ± 0.02 . 



(58) 



and the value we would expect if the network had no en- 
ergy losses is now v ~ 0.57 (recall Eqn. [32]) . The only 
noteworthy difference is that energy losses are compara- 
tively more important in this case — at about the twenty 
percent level, as opposed to ten percent in the matter 
era. There is also an enhancement of the scaling density 
that is similar to the one in the matter era. 

Fig. [5] shows the results of analogous simulations for 
slower expansion rates, namely a cosmological epoch with 
a (x t 1 / 5 and Minkowski spacetime (where the expansion 
is switched off). The former case may be cosmologically 
relevant in some toy models, for a transient epoch in the 
very early universe, but in any case it is an important nu- 
merical test for our purposes, because our analytic mod- 
eling leads us to expect that for a ex t a and a < 1/4 a 
linear scaling solution can only exist if there are energy 
losses. 

An accurate numerical evolution of this case is tech- 
nically difficult, for reasons that are well documented in 
the analogous case of cosmic strings [20] . The slower or 
non-existing expansion leads to little dissipation, which 
in turn implies that there is a large thermal bath in the 
box. In the case of strings, the optimal way to get around 
this numerical difficulty is to introduce an early period 
of gradient flow [20]. In the present case we adopted 
the simple solution of evolving the early part of these 
simulations with a matter-era expansion rate, a ex i 2 ' 3 , 
which some testing suggests is accurate enough for our 
purposes. Naturally, the slower the expansion the more 
important will be the effect of any radiation present in 
the box, and therefore the more important it is that this 
is dissipated in the matter-like evolution phase (which 
will therefore need to be longer). Our numerical tests 
suggest that for the a ex i 1 / 5 expansion rate it is suffi- 
cient to switch on the matter-like dissipation phase dur- 
ing the initial 2.5% of the simulation, but in the case of 
the Minkowski simulations this increases to about 30%. 
These values have been used for the simulations shown 
in Fig. O (which correspond the N = 3 case of the ideal 
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FIG. 6: The asymptotic values of L c /r) and u* = "/v for in- 
dividual 256 2 simulations of the ideal class of models with 
N — 3 in a slow-expansion a oc t ' 5 era (solid lines) and in 
Minkowski spacetime (dashed lines). Note that in both cases 
the early evolution has a matter-like expansion rate, in order 
to dissipate the radiation in the box which would otherwise 
contaminate the measurements. 



class of models), and the effects of the switch between 
the early and the late-time evolution are clearly visible 
in the case of the Minkowski simulation. 

Even fairly small 256 2 simulations are enough to show 
evidence for an approach to a linear scaling solution in 
both of these cases. Moreover, velocities become con- 
stant, with u* — > 1 in the Minkowski spacetime limit, 
This result therefore confirms the idea that the energy 
losses are crucially important, even for domain wall net- 
works with junctions, and that these play a fundamental 
role in the evolution of defect networks. It is possible 
that the microscopic production of coherent excitations 
of the underlying quantum field theory [2l| plays an im- 
portant role here. This issue is beyond the scope of the 
present article, but deserves further study. 



FIG. 7: The evolution of L c (top panel) and v* (bottom panel) 
as a function of the scale factor for a 1024 2 de Sitter space 
simulation of the ideal class of models with N = 10, showing 
the expected behavior: L oc a, v oc a -1 . 



C. Can the ideal model frustrate? 

Given that we are finding that linear scaling solutions 
seem to be ubiquitous for the ideal class of models, one 
might ask if this necessarily so. In other words, could 
there be some unsuspected pathological feature, hidden 
in the way the ideal class of models was designed, that 
would necessarily lead it to linear scaling and artificially 
prevent it from frustrating? However, it is easy to ex- 
clude this possibility. It is well known that domain wall 
networks should be conformally stretched (L oc a) if their 
motion is sufficiently damped [7] , for example by friction 
due to particle scattering in the early stages of their evo- 
lution or during an inflationary phase. So a simple way 
of checking that our ideal class of models can in principle 
frustrate is to show that it does reproduce this behavior. 

Performing field theory simulations of this type in de 
Sitter space involves some numerical subtleties (and, in 
particular, required some changes in our code) but can 
nevertheless be done. Fig. [7]shows the results of one such 
of 1024 2 simulation, for the N — 10 case of the ideal class 
of models. We can immediately see that the ideal class 
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FIG. 8: Comparing the scaling exponents (A) for several JV's 
in the matter era for simulation boxes of different dimensions: 
2D, 3D and 4D boxes respectively had sizes 8192 2 , 512 3 and 
128 . The error bars represent the standard deviation in an 
ensemble of 10 simulations. 



of models does behave exactly as we would expect it to 
do, namely with L ex a and v t ~ o a (HL)^ 1 ex a -1 . 
The ideal model is therefore well-behaved and it can in 
principle frustrate when this is dynamically preferred. 



D. The ideal model in other dimensions 

Let us now make a short digression to discuss the de- 
pendence of our results on the number of spacetime di- 
mensions. A simple point can illustrate why one might 
expect the possibility of such a dependence. It is well 
known [2f] that if the population probability for a given 
vacuum is larger than a critical probability p c (whose 
value depends on the lattice), that will percolate the lat- 
tice; otherwise finite vacuum bags will form. For a 3D cu- 
bic lattice p c = 0.311, so in the simplest case of a model 
with two vacua each with p — 0.5, both will percolate 
across the lattice. For a 2D square lattice p c — 0.593, so 
in the same situation neither of them percolates. Among 
other things, this means that standard domain wall net- 
works in 2D and 3D will have very different vacuum 
topologies. Nevertheless, our previous work [3, @| shows 
that this does not seem to affect the dynamical properties 
of the networks, and in particular its scaling properties. 

We now extend this result for the case of domain wall 
networks with junctions, by providing examples of the 
fact that the different dimensionality does not seem to 
affect the frustration conditions. With this aim we will 
study the behavior of some realizations of the ideal class 
of domain wall models in two and four spatial dimensions. 
(In all cases we define domain walls as the topological 
defects of co-dimension one.) 

Fig. [8] shows a comparison of the the scaling expo- 
nents we obtained in our 3D (512 3 ) matter era simula- 
tions for N = 2,8, 14, 20 (described earlier in the paper) 



with the corresponding exponents in 4D (128 4 ). A good 
agreement is found, with the slightly larger exponent and 
smaller (statistical) error bars being attributable to the 
smaller dynamic range of the 4D runs. Further evidence 
for this can be seen by performing 2D runs (8192 2 ), the 
result of which is also displayed for the case N = 2. 



E. Non-ideal models 

It is also instructive to compare the ideal mode to mod- 
els where different kinds of junctions can appear. Specifi- 
cally, we'll consider some realizations of two other classes 
of models, the Bazeia-Brito-Losano (henceforth referred 
to as BBL) and Kubotani models Kubotani [18|], Bazeia 
et al. [I9j as interesting benchmarks. In previous work 
these models allowed us to point out the key mecha- 
nisms at play, and characterize the differences between 
models with stable Y-type junctions, models with stable 
X-type junctions, and models where both types can co- 
exist. Here we will not describe these models in detail, 
but only highlight a few of their properties that will be 
relevant for comparison to the ideal class of model. We 
refer the reader to the original papers and particularly 
to Avelino et al. [llj], for a thorough characterization of 
these models. We will also revisit some of these proper- 
ties and extend our previous characterization when dis- 
cussing the role of junctions in Sect. IVI Al 

The 2-ficld BBL model has minima at the vertices of 
a square in the plane of the two fields, the orientation 
of which depends on the model parameters. The model 
allows for stable l"-type or A-type junctions depending 
on the value of the parameter which controls the tension 
of the walls connecting each pair of vacua. There are 
two classes of walls which we will denote simply by edges 
and diagonals: in the former the wall joins two neigh- 
boring minima in field space, while in the latter the wall 
joins two opposite minima. The fact that this model has 
parameters that control the ratio of energy between the 
two types of walls and consequently the type of junctions 
that are formed makes it very useful to test the energetic 
and geometric considerations discussed in Sect. Ill Dl 

For models with three scalar fields there is an analo- 
gous BBL model but an alternative is provided by the 
(perturbed 0(3)) Kubotani model. Again there are two 
branches for the minima, which will be placed either at 
at the vertices or at the centers of the faces of a cube in 
the space of the scalar fields. In the former case there are 
three kinds of walls, which for obvious reasons we can re- 
fer to as edges, external diagonals and internal diagonals 
while in the latter there are two kinds of walls, which we 
can refer to as edges and axes. Again the choice of the 
model parameters determines what type of junctions will 
be present: in addition to the parameter ranges where 
only X-junctions or Y-junctions are stable, now there is 
also a range of parameter where both are stable and co- 
exist. The 3-field BBL and Kubotani models are quite 
similar, and it is possible to exhibit a correspondence be- 
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FIG. 9: The values of L c /r] as a function r/ for the ideal class 
of models (solid lines) with N ranging from 2 to 20 (from top 
to bottom respectively) and for the BBL model with 2 and 
3 fields (dashed top and bottom lines respectively). Details 
of these models can be found in [ill ]. For both BBL models 
the simulations have been with the model parameter choice 
er 2 — —3/10. All lines correspond to an average of ten 128 3 
simulations. 



tween both potentials, although this is not one-to-one. 
This implies the Kubotani-model simulations can also be 
interpreted as BBL model simulations, although the con- 
verse need not be true. 

Our present intention is to compare the ideal model 
with these other models from the point of view of pro- 
ducing frustrated networks. We have performed series 
of 128 3 matter-era simulations of the BBL model in the 
case where only A-type junctions arc allowed. In both 
the 2-field and the 3-field case for the parameter choice 
er 2 = —3/10 will lead to this scenario — again, refer to 
[111 ] for a discussion of this parameter. 

In Fig. [5] we compare the values of L c jr\ as a function 
r\ for the ideal class of models, with N ranging form 2 
to 20, with the aforementioned BBL models with stable 
X-type junctions. As expected, we find that the 3-field 
BBL model gives a better result than the 2-field one, but 
there is no significant improvements with respect to the 
ideal class of models. The no frustration conjecture is 
evident for all models. 



V. FURTHER MODELING 

In this section we will provide an analysis of some of the 
above results based on the simple analytic model which 
we introduced in Sect. Ill CI The model was derived and 
extensively described in [7|, a nd is analogous to an earlier 
model for cosmic strings [22J . There it was also shown 
that it provides a good description of the evolution of the 
large-scale properties of standard domain wall networks. 
We also discuss the behavior of this type of networks in 
the ideal limit N — > oo. 



A. The low-A case 

For standard domain walls a reasonable fit to the nu- 
merical simulations, both in the radiation and in the mat- 
ter eras, is obtained for the following choice of model 
parameters 



1 



co ~ k /2 . 



(59) 



Note that since this is a fairly simple phenomenological 
model we won't include any error bars in the model pa- 
rameters; these would typically be at the ten to twenty 
percent level. 

Let us start by concentrating on the simulations of the 
case N = 2. Since the late-time evolution is found to be 
close to the linear scaling regime, we can approximately 
use the scaling relations given by Eqns. [551 USB an d we 
then find that the analytic model still provides a good fit 
to the data, although the phenomenological parameters 



fc 2 ~0.7, c 2 ~fc 2 /2 



(60) 



This holds for both the matter and radiation eras, al- 
though the fit is better in the former case than in the 
latter. This bias towards the matter era simulations in 
the choice of fit parameters is deliberate, since we have 
done larger simulations in the matter era, and the faster 
expansion rate (relative to that the radiation era) also 
means that these simulations should be less affected by 
the thermal bath in the box. 

As was to be expected, we find parameter values that 
are significantly smaller than in the case of standard do- 
main walls. The presence of junctions means that the 
walls tend to be straighter, and hence k (which is to 
some extent has the phenomenological role of a curvature 
parameter) must be smaller than in the standard case. 
Nevertheless, note that it is still quite close to unity — 
our earlier work |1(J, QjJ shows that it would need to be 
several orders of magnitude smaller for frustration to be 
viable. 

Now let us turn our attention to the behavior of the 
scaling properties as we increase the number of fields. We 
have already pointed out when presenting the simulation 
results that the fact that the scaling velocities are inde- 
pendent of the number of fields is not really surprising. 
The key features of the ideal class of models are that all 
the minima of the potential are equally spaced and all 
the domain walls have the same tension (a model with 
N fields will have N + 1 minima). These in turn nec- 
essarily imply [101 ] that the junctions will always be of 
Y type, regardless of the number of fields. Naively, on 
therefore expects that in a statistical sense all walls of a 
given network (in a given cosmological epoch and with a 
given number of scalar fields) probe a similar region of 
the potential, and will have the same basic properties, 
such as the characteristic scale and root-mean squared 
velocity. In these circumstances, recalling our discussion 
in Sect. Ill Al and provided we can assume that the junc- 
tions themselves do not have any crucial dynamical role, 
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then the local dynamics of a single domain wall will ef- 
fectively be roughly the same regardless of the number of 
fields. We will return to the issue of the dynamical role 
of the junctions in the following section. 

As an aside, recalling that the equation of state of a 
domain wall gas is given by 



P 2 

w=~=-- 

P 3 



(61) 



our results on the network velocities allow us to predict 
the following equations of state for scaling domain wall 
networks with junctions in the radiation and matter eras 
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1-Urad r ^ 1 i r ' tVmat ^ „ . ■ 



(62) 



so the correct value is closer to w ~ —1/2 than to the 
naively expected w = —2/3. We emphasize that for the 
(ideal) class of models we are considering, these will hold 
for any number of fields. It's also clear that these are in 
conflict with the existing observational bounds. 



B. TV dependence and the ideal limit 

Changing the number of fields N will change the net- 
work's characteristic length, and hence its density. Look- 
ing at the plots for the scaling density parameter e = L/t 
in the matter and radiation eras, it is strikingly clear that 



e low {N) ex 



N 



(63) 



provides a good fit to the data. If this fit were to hold 
for any number of fields then it would obviously imply 
that e^v — > as N — > oo. In other words, even for a fixed 
defect tension one could freely tune the energy density of 
the network at a given subsequent epoch to take any value 
we wanted, simply by choosing the appropriate number 
of fields — although this would be a clear form of fine- 
tuning. However, it is clear that Eqn. [53] is not valid for 
all values of N, and specifically in the limit N — > oo. In 
fact in the large N limit we expect a behavior of the form 



Zhigh(N) =A + 



B 

N 



(64) 



corresponding to a constant asymptotic density. The rea- 
son is clear. On one hand, L/(vi) has to be of order unity 
and we have seen that the characteristic velocity of the 
domain walls, v, is rather large and appears to be inde- 
pendent of N. On the other hand the probability that 
the collapse of a given domain results in the fusion of ad- 
jacent domains is proportional to 1/iV, for large N. The 
fusion of adjacent domains will change the domain wall 
density with respect to that of the N — + oo limit by a 
factor proportional to the probability of such events, if 
N is large, which leads to Eqn. [53] in this limit. 

However, this /V-dependence overestimates e (thus un- 
derestimating the density) for small numbers of fields, 



but an adequate fit can be obtained for the larger values 
we simulated. Specifically, in the range 10 < N < 20 we 
find 



B n 



1.0, B r 



1.2, 



(65) 



respectively in the matter and radiation eras, while in 
both cases 



A -0.06, 



(66) 



which would therefore characterize the asymptotic scal- 
ing law. 

Regardless of the detailed dependence on the number 
of fields, for sufficiently large energy densities the wall 
network will at some point become the dominant energy 
component, and the subsequent evolution of the universe 
will be dramatically altered — see Q for a brief discussion 
of this scenario. And even before this stage is reached, 
a sub-dominant but significant contribution to the total 
energy density might have consequences that will conflict 
with existing observations — we will return to this point 
in Sect. IVII1 Before that, however, we will consider what 
changes (if any) are to be expected in the above results 
if the junctions are dynamically important. 



VI. THE ROLE OF JUNCTIONS 

The types and properties of junctions that can form 
will depend both on energy considerations and on the 
topology of the potential's minima in field space. Since at 
least n real scalar fields are required to produce a defect 
of co-dimension n, it follows that a junction of the respec- 
tive type will correspond to an n-dimensional field space 
configuration, connecting at least [n + 1) minima. Note 
that one can have configurations where defects intersect 
but there are no physical junctions (a trivial example is 
the case of decoupled fields). One can then discuss ex- 
amples of particular models, as we will do below. We will 
then briefly address the issue of the role of the junctions 
themselves in the dynamics of the defect networks. 



A. Examples of junction hierarchies 

In what follows we will discuss the possible junction 
hierarchies in three classes of models we have been con- 
sidering: in addition to the ideal class we have been dis- 
cussing for most of this article, we will again use the BBL 
and Kubotani models [la 13 as interesting benchmarks. 
The present brief discussion can be seen as an extension 
of that in [ll| (we leave a more detailed analysis, par- 
ticularly of energetic arguments, for future work). Here 
as in that article, the point we wish to emphasize is that 
despite very simple rules for the formation of junctions, 
a combination of topological, geometric and energetic ar- 
guments can lead to quite different behaviors, depending 
on the details of the particular model. As one might 
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expect, the ideal class of models will turn out to have 
the simplest and physically clearer phenomenology. No- 
tice that for the sake of clarity the discussion will focus 
on the case of 3 spatial dimensions — however, it's simple 
enough to repeat for other dimensions. 

In the 2-ficld BBL model there are always 4 min- 
ima, placed at the vertices of a square. Domain walls 
form, and intersect at string-type junctions. Depending 
on the model parameters, energetic configurations will 
imply that either X-type or Y-type wall junctions will 
be stable (that is, with either 4 or 3 walls meeting at 
each junction) — see [ll[ for a discussion of both cases. 
There is no situation in which both are simultaneously 
stable. Notice that in the former case there is only one 
type of string that can be formed, while in the latter 
there are 4 types, although they all have the same en- 
ergy. Also notice that in this model there are physically 
no monopole-type junctions — which is obvious given that 
there are only two scalar fields. 

In the positive branch of the Kubotani model (which 
corresponds to the intermediate range of the BBL 3-field 
model) , there are 3 scalar fields and 8 minima which oc- 
cupy the vertices of a cube. Walls form string-type junc- 
tions which, by energetic arguments are always of X-type 
(the Y- types being unstable in this regime), so they form 
along the edges of the cube. Hence 4 walls meet at a 
string-type junctions, but now the hierarchy continues 
one further step, and 6 strings meet at a monopole-type 
junction. Notice that there is only one type of (stable) 
monopole, where 12 walls effectively meet (and not 24 — 
each wall is effectively contributing to two strings). This 
is illustrated in the top panel of Fig. [lOj 

The negative branch of the Kubotani model (which 
is also the negative range of the BBL 3-field model) is 
slightly more subtle. Here there are 6 minima at the 
vertices of an octahedron. Energetic arguments show 
that we now have walls with Y-type junctions form- 
ing along the octahedron's edges, so 3 walls can meet 
at a string-type junctions and 8 strings can meet at a 
monopole-type junction. This is illustrated in the bot- 
tom panel of Fig. [TQl However, in this case there are 
also other types of junctions allowed, both for strings 
and monopoles. It turns out that energetic considera- 
tions also allow for walls to form X-type junctions, as 
well as Y-type junctions involving the axis sector — both 
of these have in common the fact that they are con- 
figurations in which one of the fields vanishes. From 
this it follows that there are also monopole-type junc- 
tions formed at the intersection of 5 strings (these have 
4 strings formed from edge sector Y-type junctions and 
one formed from an axis sector X-type junction) or at the 
intersection of 4 strings (this involves 2 strings formed 
from edge sector Y-type junctions and 2 strings formed 
from axis sector Y-type junctions). Pictorially these 
monopoles correspond to topological configurations in- 
volving one half or one quarter of the octahedron respec- 
tively (or equivalently, square-based and triangular-based 
pyramids), while the 8-string monopole corresponds to 




FIG. 10: Two possible junction hierarchies for the Kubotani 
model in three spatial dimensions. In the so-called positive 
branch (top panel) 4 walls meet at a string-type junctions and 
6 strings meet at a monopole-type junction. In the negative 
branch (bottom panel) 3 walls can meet at a string-type junc- 
tions and 8 strings can meet at a monopole-type junction. In 
the latter case other types of junctions are also aallowed — 
refer to the main text for a discussion. 



the full octahedron. The relation between the energies of 
the different monopole-type configurations is not imme- 
diately clear — in particular, it is conceivable that some 
are unstable and will tend to decay. 

The positive range of the BBL 3-field model (which has 
no correspondence in the Kubotani model) is even more 
complicated. It has 8 minima at the vertices of a cube 
(just like the intermediate range) but now energetic argu- 
ments imply that the walls meet at Y-type junctions (the 
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FIG. 11: The junction hierarchy in the ideal class of models, 
in three spatial dimensions. String-type junctions form at the 
intersection of 3 walls, and monopole-type junctions form at 
the intersection of 4 strings. 



X-type ones are unstable in this parameter range). Hence 
string-type junctions form at the intersection of 3 walls, 
but one can have three physically distinct classes of junc- 
tions: they can involve 3 walls from external diagonals, or 
2 walls from edge sectors and one from an external diag- 
onal, or one wall from each of the sectors (edge, external 
diagonal and internal diagonal). We can refer to each of 
these strings as Type I, II and III respectively, according 
to the number of different wall types involved. Assuming 
that all three classes are allowed, monopole-type junc- 
tions can form at the intersection of 4 strings, and again 
there will be 3 possible types of monopoles. Interestingly, 
they all must necessarily involve a combination of differ- 
ent types of the wall junctions: one can have monopoles 
at the intersection of three Type II and one Type I string, 
or two Type II and two Type III strings, or finally one 
Type I, one Type II and two Type III strings. In any 
case, 6 walls effectively meet at a monopole, as opposed 
to the 12 of the intermediate range discussed above. Of 
these 6, there will be at least 2 walls from edge sectors 
and 2 from external diagonals. The remaining 2 will be 
of different types: in the first case we discussed we have 
one edge and one external diagonal, in the second case we 
have one edge and one internal diagonal, and in the third 
case one external diagonal and one internal diagonal. It 
would be interesting to study the relative energies of all 
these configurations in more detail, and hence their sta- 
bility, both analytically and numerically. We expect that 
the three monopole cases we listed abovr are ordered by 
increasing energy. 

Finally, in the ideal model energetic arguments im- 
ply that walls always form Y-type junctions, and with 
p scalar fields there will be (p+ 1) minima evenly spaced 
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FIG. 12: The averaged velocities (with statistical error bars) 
of 1024 3 matter simulations of the ideal class of models with 
N — 2 (dashed lines) and N = 3 (solid lines). The velocities 
are plotted as a function of the simulation time, and in each 
case the averages are over 10 simulations. 



on the surface of a p-dimensional sphere, from which we 
will produce defects of co-dimension up to p. Hence in 
three spatial dimensions string- type junctions form at the 
intersection of 3 walls (provided there are at least 2 fields) 
and monopole-type junctions form at the intersection of 
4 strings (provided there are at least 3 fields). An illus- 
tration of this configuration is in Fig. [TTJ In models with 
additional dimensions there can be a much deeper hier- 
archy: the rule is that co-dimension (p + 1) defects will 
form at the intersection of (p + 2) defects of co-dimension 
p (for p > 0) . 



B. Junctions and network dynamics 

From the above discussion of the ideal model it's note- 
worthy that the case N = 2 of the ideal model is special 
in that it contains no true monopoles, while those with 
N > 2 do have them. (Trivially, the case N = 1 has 
neither monopoles nor strings, only plain domain walls.) 
It is therefore interesting to compare this case with say 
N = 3: if the differences are larger than expected from 
the discussion in the previous sections, this could be an 
indication of an effect of the monopoles on the dynamics 
of the network. In order to have a closer look at this, we 
have carried out 1024 3 matter-era numerical simulations 
of the N = 2 and N = 3 cases. 

Fig. [T2] shows the evolution of the velocities in the two 
cases (averaged for series of 10 runs of each type). The 
fact that the velocities are statistically the same (notice 
the very narrow range of the y axis) naively suggests that 
the dynamical effect of the monopoles is negligible. Nev- 
ertheless, it is of interest to consider what changes would 
be expected if the junctions are dynamically important. 
It's clear that this will increase the effective equation of 
state of the defect network and therefore it will not be 
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any help for dark energy and frustration scenarios, but 
it might have other interesting consequences. 

Very simple models can be used to study this issue — 
even 2-field models and 2D simulations are sufficient to 
phenomenologically study the effect of the junctions. The 
idea is to consider walls which have fixed tension and 
thickness, but also a junction energy that is a function of 
a single tunable parameter. One way to implement this 
would be a perturbed version of the N = 2 ideal model. 
An alternative is a to take a perturbed Z 3 model, namely 



L = \\y>-\){<? 



1) 



6{v 2 ~l) 2 



(67) 



where 9 is a positive parameter which controls the junc- 
tion energy. However, notice that in order to keep a fixed 
thickness we must modify this to 



L = ±\*f(0)(<p 3 



i)(^-i) + e^ 2 -iy 



(68) 



where f(9) is a function to be determined numerically. 
One can then study the scaling properties and/or the 
effective equation of state as a function of the parameter 
9. This will be left for future work. 



VII. DISCUSSION AND CONCLUSIONS 

In this paper we have built upon earlier work [lOJ, I 111 . 
Il6| and presented compelling evidence that domain wall 
networks cannot be the dark energy. Note that in or- 
der to be able to rule out the domain wall scenario for 
dark energy a very large and rich class of models had to 
be analyzed in detail. This led us to develop a model 
best suited for frustration (the 'ideal' model). We have 
shown that even this model fails to produce a frustrated 
domain wall network. Current observational constraints 
using cosmic microwave background and supernova data 
already strongly disfavor w = —2/3 as the equation of 
state of a single dark energy component [5|. However, 
we should bear in mind that these observational results 
are dependent on strong priors. Even if we take them for 
granted and accept that domain walls alone cannot be 
the dark energy, die-hards might argue that they could 
still make a significant partial contribution. Our results, 
however, exclude even that rather more contrived possi- 
bility. 

More generally, we have also provided a quantitative 
characterization of the cosmological evolution of the ideal 
class of domain wall networks with junctions. These have 
the interesting feature that their characteristic velocity 
appears to be independent of the number of scalar fields 
N. We have also shown how a simple analytic toy model 
can provide an adequate description of these models and 
discussed how the scaling density depends on the num- 
ber of fields, both for the numerically simpler case of a 
small number of fields and for the physically more inter- 
esting case N — *■ oo. We have also discussed the possible 



hierarchies of junctions in the ideal class of models, con- 
trasted it with what is found in other models and briefly 
addressed the role of the junctions in the network dy- 
namics. 

We are now in a position re-examine the possible cos- 
mological roles of domain wall networks and the corre- 
sponding bounds on relevant model parameters in the 
light of our findings. The Zel'dovich bound |3J] tell us 
that the energy scale T$b of the symmetry breaking 
phase transition that produced the domain walls has 
to be smaller than about 1 MeV. However, the clas- 
sic derivation of the Zel'dovich bound implicitly assumes 
the linear scaling solution with roughly one defect per 
Hubble volume. As we saw above this is not in general 
a valid assumption, and in particular it is quite wrong if 
the domain walls are to have any chance of providing an 
interesting contribution to the dark energy. 

We have shown that L < lOkpc at the present time 
and consequently LH < 1CP 6 if the domain walls are to 
provide a significant contribution to the dark energy. Of 
course since the characteristic scale of the domain walls 
has to be smaller than the Hubble scale by a factor of at 
least 10 -6 then the domain wall tension also has to be 
smaller by a similar factor (with respect to the Zel'dovich 
bound) in order that domain wall density does exceed the 
critical density. This translates into a bound on the scale 
of the symmetry breaking phase transition that produced 
the domain walls of 



T SB < 10 keV . 



(69) 



Finally, let us point out that our massively parallel do- 
main walls code, which has been gradually optimized and 
generalized to be able to deal efficiently with multi-field 
scenarios, can have other uses beyond the study of do- 
main wall networks (with or without junctions). With 
suitable changes it can be used to phenomenologically 
study cosmic (super)strings, which are generically ex- 
pected to have junctions, and also to look into a num- 
ber of string landscape scenarios, where the key feature 
is the presence of a large number of (usually) coupled 
fields. Some of these topics will be the subject of future 
publications. 
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